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Abstract 

The exact solution for the motion of a test particle in a non-spherical 
polar orbit around a Kerr black hole is derived. Exact novel expressions 
for frame dragging (Lense-Thirring effect), periapsis advance and the or- 
bital period are produced. The resulting formulae, are expressed in terms 
of Appell's first hypergeometric function F\, Jacobi's amplitude function, 
and Appell's -Fi and Gaufi hypergeometric function respectively. The ex- 
act expression for frame dragging is applied for the calculation of the 
Lense-Thirring effect for the orbits of S-stars in the central arcsecond of 
our Galaxy assuming that the galactic centre is a Kerr black hole, for 
various values of the Kerr parameter including those supported by recent 
observations. In addition, we apply our solutions for the calculation of 
frame dragging and periapsis advance for stellar non-spherical polar or- 
bits in regions of strong gravitational field close to the event horizon of 
the galactic black hole, e.g. for orbits in the central milliarcsecond of 
our galaxy. Such orbits are the target of the GRAVITY experiment. We 
provide examples with orbital periods in the range of lOOmin - 54 days. 
Detection of such stellar orbits will allow the possibility of measuring the 
relativistic effect of periapsis advance with high precision at the strong 
field realm of general relativity. Further, an exact closed form formula for 
the orbital period of a test particle in a non-circular equatorial motion 
around a Kerr black hole is produced. We also derive exact expressions 
for the periapsis advance and the orbital period for a test particle in a 
non-circular equatorial motion in the Kerr field in the presence of the cos- 
mological constant in terms of Lauricella's fourth hypergeometric function 
Fd. 
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1 Introduction 
1.1 Motivation 



Most of the celestial bodies deviate very little from spherical symmetry, and 
the Schwarzschild spacetime is an appropriate approximation for their gravi- 
tational field [2]. However, for some astrophysical bodies the rotation of the 
mass distribution cannot be neglected. A more general spacetime solution of 
the gravitational field equations should take this property into account. In this 
respect, the Kerr solution [3] represents, the curved spacetime geometry sur- 
rounding a rotating mass [4]. Moreover, the above solution is also important 
for probing the strong gravitational field regime of general relativity [5] . This 
is significant, since general relativity has triumphed in large-scale cosmology 
[8, 7, 9, 10], and in predicting solar system effects on planetary orbits like the 
perihelion precession of Mercury with a very high precision [1, 6] ^. 

Due to their significance for black hole physics and for describing orbits 
around a rotating mass, exact solutions of Kerr geodesies with and without 
the contribution from the cosmological constant have been derived and studied 
extensively in recent years [11, 12]. However, few important classes of possible 
orbits were left out of the discussion in [11, 12]. 

An important such class is the case of timelike non — spherical polar Kerr 
orbits. An essential property of such non-equatorial orbits in Kerr spacetime 
is the Lense-Thirring phenomenon or frame dragging: the Kerr rotation adds 
longitudinal dragging to the relativistic orbital precession in the Schwarzschild 
spacetime. Thus the orbital plane for such geodesies in Kerr spacetime is not 
fixed but rather is dragged by the rotation of central mass in sharp contrast with 
orbital solutions in the static Schwarzschild case. It is the purpose of the present 
work to derive the first exact solution of this important class of orbits in closed 
analytic form. Our exact analytic solutions provide closed form expressions for 
the orbit, Lense-Thirring effect, periapsis advance and the corresponding orbital 
period. Further our exact results allow to precisely calculate the Lense-Thirring 
period (frequency) for these non-zero eccentricity orbits. By construction our 
solutions are valid at all regions of gravitational strength from the strong to 
the weak gravitational field regime. Thus they are particularly important for 
probing the strong field realm of general relativity. 

Another important class studied in this work for the first time is the case 
of equatorial non — circular timelike orbits in the presence of the cosmologi- 
cal constant A. In particular, we derive the exact solutions for the periapsis 
advance and the orbital period for a test particle in a non-circular equatorial 
motion in the Kerr-de Sitter gravitational field. The derived exact formulae in 
closed analytic form are important for the precise determination of the com- 
bined effects of A and the rotation of central mass on planetary orbits as well 

^See also about the BepiColombo science mission on Mercury 
http://sci.esa.int/home/bepicolombo/. 



2 



as for investigating the effect of the cosmological constant in black hole physics 
e.g. accretion physics for non-zero eccentricity trajectories. 

Motivated by the observational evidence that the galactic centre region Sgr 
A* , of Milky Way. is the best candidate for a superniassive rotating black hole 
[15, 16], we applied our solutions for non-spherical polar Kerr orbits in two 
situations a) for test particle stellar orbits near the outer horizon of the galactic 
centre black hole (strong field regime) b) for calculating the frame dragging of 
the observed orbits of S-stars in the central arcsecond of our galaxy (weak field 
regime). 

Indeed observations of the galactic nucleus at the near-infrared and X-rays 
have provided a quite compelling evidence for the existence of a supermassive 
black hole of mass Mbh = (3 — 4) x IO^Mq at the galactic centre, located at 
a distance of approximately 8Kpc [13, 14, 15, 16, 18] {Mq denotes the mass of 
the Sun ). 

Further observations of near-infrared periodic flares have revealed that the 

central black hole is rotating so that the spacetime region near Sgr A* is de- 
scribed by the Kerr geometry rather that the Schwarzschild static geometry, 
with the Kerr parameter [15] 

= 0.52 (±0.1, ±0.08, ±0.08) (1) 



GMbu/c 

where the reported high-resolution infrared observations of Sgr A* revealed 'qui- 
escent' emission and several flares. These flare events last for about an 1 h, and 
their light-curves show significant variations on a typical scale of 17min. In 
the above equation J ^ denotes the angular momentum of the black hole (The 
error estimates here the uncertainties in the period, black hole mass and dis- 
tance to the galactic centre, respectively; G is the gravitational constant and 
c the velocity of light.) This is half the maximum value for a Kerr black hole 
[24] . Observation of X-ray flares conflrmed that the spin of the supermassive 
black hole is indeed substantial and values of the Kerr parameter as high as 
oiGaiactic = 0.9939 have been recently reported [16, 27]. Let us also mention that 
a hump of the locally non-rotating frames (LNRF) velocity profile has been 
found for accretion discs orbiting rapidly rotating Kerr black holes with spin 
a > 0.9953 for Keplerian discs [27] and a > 0.99979 for marginally stable discs 
[28] 3. 

The material of this paper is organized as follows. Section 2 and 3 deals 
with polar geodesies in a Kerr background while section 4 discusses equatorial 

geodesies in a (cosmological) Kerr geometry. In section 2 we derive the first ex- 
act orbital solution that describes the motion of a test particle in a non-spherical 



^ J = ca where a is the Kerr parameter. The interpretation of ca as the angular momentum 
per unit mass was first given by Boyer and Price [25] . In fact, by comparing with the Lense- 

Thirring calculations [26] they determined the Kerr parameter to be; a = — where O 

and I denote the angular velocity and radius of the rotating sphere. 

^Such humpy equatorial orbital velocity profiles for stable, circular orbits, occur in a small 
radial range close to but above the marginally stable equatorial circular geodesic of the black 
hole [27, 28). 
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polar orbit in the Kerr gravitational field. More specifically, the exact expression 
in closed analytic form for the amount of frame dragging (Lense-Thirring) effect 
that a test particle undergoes is given in Eq.(12) of subsection 2.0.1, and its rel- 
ativistic periapsis advance is given in Eq.(29) of subsection 2.0.2. In subsection 
2.1 we derive an exact expression for the orbital period of a polar non-spherical 
Kerr orbit Eq.(33) which is important for defining the Lense-Thirring period 
Eq.(37). We subsequently apply our results to an important system for which 
observational data is available. Namely i) in 2.2 we calculate the frame dragging 
for the observed orbits of S-stars in the central arcsecond of our Galaxy assum- 
ing that the galactic centre is a rotating Kerr black hole. We calculated the 
Lense-Thirring periods for various values of the Kerr parameter and compared 
the result with estimated ages of the S-stars from spectroscopic studies, ii) We 
calculate in section 3.0.1 the frame dragging and periapsis advance for polar 
non-spherical Kerr orbits in the central milliarcsecond of our Galaxy. In section 
4 we derive the exact solution for the relativistic periastron advance, Eq.(53), 
for the equatorial orbital non-circular motion of a test particle in the Kerr field 
in the presence of the cosmological constant. We also find the orbital period, 
Eq.(58), for equatorial non-circular orbits in the asymptotically flat Kerr black 
hole. 

The angular resolutions required for measuring stellar orbits in the central 
milliarcsecond of Milky Way will be achievable by the General Relativity Anal- 
ysis via Vlt InTerferometrY (GRAVITY) experiment using Very Long Baseline 
Interferometry (VLBI) techniques [17]. The detection of such stellar cusp orbits, 
predicted by the theory will allow a precise measurement of the relativistic 
effect of periapsis advance since the orbital periods are small and therefore will 
allow a test of the precise theory developed in this work, based on the exact so- 
lutions of the equations of General Relativity at the strong-field realm. A com- 
pletely unexplored region of paramount importance for fundamental physics. 
We present examples of non-zero eccentricity trajectories with orbital periods 
in the range lOOmin-54 days. In a particular solution that will be presented 
in section 3.0.1 with orbital period of 53.6 days the phenomenon of relativistic 
periapsis advance is 19-8^, already a very large effect, which therefore should 
be easily measurable with a high precision, after say one year of observation. 

In a scries of appendices we present some of our formal calculations as well 
as some background on Lauricella's hypergeometric function, the differential 
equations it obeys as well as its integral representation. 

Taking into account the contribution from the cosmological constant, the 
generalization of the Kerr solution is described by the Kerr -de Sitter metric 
element which in Boyer-Lindquist (BL) coordinates ^ is given by [30, 31]: 

^Stellar counts expect an apparent magnitude in the range 17.5 < rriK < 19-5 for stars 
within this region of the galactic centre. 

^ These coordinates have the advantage that reduce to the Schwarzschild solution with a 
cosmological constant in the limit a —* 0, see [29]. 
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where 



ds"^ = (cdt - asin^ 0d(j)f - ^dr^ - ^dO'^ 



Ae sin^ 9 



{acdt - {r^ + a^)d(l>) (2) 



^ 2\f 2 2\ 2GA'/bH 



A, := {l--r')(r'+a') 

Ae := 1 + — cos^ e 



a 



2 2,2 2 /1 

_ .— X -I — , p := r + a COS 

o 

(3) 

and A denotes the cosmological constant, a is the Kerr (rotation) parameter 
and Mbh is the mass of the black hole. 

The timelike geodesic equations in the presence of the cosmological constant 
have been derived in [11] by solving the Hamilton- Jacobi partial differential 
equation with the method of separation of variables: 



p^^ = =^ (a^sin^ e-L) + ^ \{r^ + o?)E - ah] 

^ ds Ae sin^ 9^ ' A^ ' ^ 

2 dt _ S2(r2 + a2) [(r^ + a^)E - ah] aE^{aE sin^ 9 - L) 



ds Ar A, 

dr 
ds 



p'^ = ±Vr 



p':r. = ±v^ (4) 



where 



R := E^[{r^+a^)E-aLf -Ar{r^ + Q + E\L-aEf) 

e := \Q + {L- aEfa? - cos' 9] Ae - ~? ("^^^^"'^ " -^)' (5) 

sin 6 

The constants of integration E, L arc associated with the isometries of the 
Kerr metric. Carter's constant of integration is denoted by Q. 

The geodesic differential equations (4), (5) reduce to the ones found by 
Carter [34] when we set the cosmological constant to zero. 
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2 Exact solution of non-spherical polar timelike 
Kerr geodesies 

2.0.1 Frame dragging effect for polar non-spherical bound orbits 

In this section we shall derive the exact solution for the important class of 
non-sphcrical polar orbits of a test particle in the gravitational field of a Kerr 
black hole. We shall first determine the change in the azimuthal angle (f) after a 
complete radial oscillation for the case of non-spherical polar timelike geodesies, 
assuming a vanishing cosmological constant. This will allow us to calculate the 
resulting frame-dragging (Lcnsc-Thirring effect). Thus we generalize our earlier 
results in which the exact solutions of the spherical polar timelike geodesies 
have been derived and applied for the determination of frame dragging (Lense- 
Thirring) effect for stellar orbirs close to the galactic centre, and for satellite 
orbits around the Earth [11, 12]. It has been shown [40] that a necessary con- 
dition for an orbit to be polar (meaning to intersect the symmetry axis of the 
Kerr gravitational field) is the vanishing of the parameter L, i.e . L = ®. The 
relevant differential equation for the calculation of frame dragging is 

dr aVR ^ ' 

where the quartic polynomial R, for non-spherical polar Kerr orbits, is given 

by (5) with A= 0,L = 0. Also A := + - 2GMBHr. As in [12] we will 
use partial fractions for integrating equation (6) from periapsis distance rp to 
apoapsis distance ta- Thus we have 

# _ ^+ + ^- 



dr {r-r+)\/R {r - r-)\/R 



+ 



(r - r+)^{E^ - l)(r - a){r - /3)(r - ^){r - S) 

AP 

(r - r-)^/{E^ - l)(r - a)(r - /3)(r - 7)(r - 6) 



(7) 



We denote the real roots of the quartic polynomial by a, /3, 7, 5, a > /3 > 7 > (5 
and organize all roots in ascending order of magnitude as follows 

(8) 



®Such orbits are also characterized by a non- vanishing value for the hidden first integral 
represented by Carter's constant Q. 
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where ap = a^, = a^+i, — ap,+2-, ca = a^-i, i = 1,2,3 and we have that 
Q^/i-i > C(fi-2 > ctn-s- By applying the transformation ^ 

z=i'-^^^ (9) 
w r' - 

where 

a;=^^i^^ (10) 
an - cin+2 

and after a dimensionlcss variable r' through r = r'^^Mp^ has been introduced, 
we can bring the radial integral which determines the amount of frame dragging 
that a non-spherical polar Kerr orbit undergoes in the general theory of relativity 
(GTR) 



/ r- (11) 



irp 

into the familiar integral representation of Appell's first hypergeometric function 

of two variables Fi . 



Indeed we obtain 



H+ ^2' ' 2' ' ; r(2) 

^3/2^p /3 1 . 2\ r(| r(i) 

0.1/2 /I 1 „ 2\ r^i) 



(12) 



The function Fi{a, (3, /?', 7, .t, y) is the first of the four Appell's hypergeometric 
functions of two variables x,y [38] ^, 



00 C30 



which admits the following integral representation 

£ u'^-\l-uy-'^-\l-ux)-^{l-uy)-^'du = rW-T- cx) ^^^^^ ^, ^ 

^ (14) 



^Wc have the correspondence: = a, a^+i = I3,a^^2 = 7iQn-i = ^^-2 = 'f''±,o.ij,—'i 

S. 

*The expression (A, k) = A(A + 1) • • • (A + /« — 1), and the symbol (A, 0) represents 1. 
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The double scries converges when \x\ < 1 and \y\ < 1. The above Euler 
integral representation is valid for Re(Q:) > and Re(7 — a) > 0. Also 
T{p) = xP~^e~^dx denotes the gamma function. In addition 

= V(l - E^){P - r'±)V^^VJ^ (15) 

and 

The radii of the event horizon and the inner or Cauchy horizon, r+ , r_ respec- 
tively of the black hole are located at 



r± = ^±^[-^) (17) 

In equation (16) rj_ := oli^n i o, are the dimensionless horizon radii and spin of 

the black hole; respectively. 

The moduli, introduced in (12), are given in terms of the roots of the quartic 
and the horizon radii of the rotating black hole by the expressions 
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o-ti - Q-n+i oin-i - an+2 ^ g - /3 (5 - 7 



(18) 



Equation (12) represents the first exact formula of frame dragging (Lense- 
Thirring effect) that a test particle in a non-spherical polar bound orbit around 
a Kerr black hole undergoes. It represents a generalization of our previous 
results for test particles and photons in spherical polar orbits [11, 12]. 



2.0.2 Periapsis advance for non-spherical polar orbits 

In this subsection we shall investigate the periastron advance for a non-spherical 
bound Kerr polar orbit, assuming a vanishing cosmological constant. We shall 
first obtain the exact solution for the orbit and then derive an exact closed form 
formula for the periapsis precession. The relevant differential equation is 

where in (19), Q and R are quantities defined, respectively, in (5) with A= 
0, i = 0. Now applying the transformation (9) on the left hand side of equation 
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(19) we get 

dr 1 f dr' 



1 f dzy/uj 



GMbh 



where 



V(l - E^)^{a - - 6)^z{l - z){l - kH) 
a — P S — J 5 — 7 



CO- 



(20) 



(21) 



a — 7(5 — /? 5 — (3 
By setting, z = x'^, we obtain the equation 

[ dx _ Vl - E'^y/'^r^y/^^ r d0 

J ^/{i-x^){i- K^x^) " i 7i 

Using the idea of invertion for the orbital eUiptic integral on the left hand side 
we obtain 

^^^(^^EEvSlLZJsM^J^,^ (23, 
In terms of the original variables we derive the equation 



^ 1^ 2 j ve' ^ ; GMbh 

a-7 I 2 J y^' 



(24) 



Equation (24) represents the first exact solution that describes the motion of a 
test particle in a polar non-spherical bound orbit in the Kerr field in terms of Ja- 
cobi's sinus amplitudinous elliptic fimction. The fimction sn^(?/, k^) has period 
2K{k^) = ttF (i, 1, k^) which is also the period of r The Gaufi hypergeo- 

metric function F(ai,/3i, 7]^, a;) is equal to 1 + "^'J^^^ x + ^"2"*^^ fa+i)^"^^ 

[41]. This means that after one complete revolution the angular integration has 
to satisfy the equation 

de I f 95- 4 1 1 TT /I 1 . 



V& VQJ Vl-«'2sin2* VT^^V^^VP^^ V2'2 

(25) 

where the latitude variable 

*:=7r/2-6l (26) 



^The roots of the quartic axe organized as: a^i > a^+i > a;^+2 > a^^-a, and we have the 
correspondence a = afj,,l3 = 0:^+1,7 = = 

i^The roots of the quartic and the angular integration are understood to be dimensionless 
in this equation. 

^^K{k'^) denotes the complete elliptic integral of the first kind. 
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and 



has been introduced. Equation (25) can be rewritten as 

dx' 

"7(^-x'2)(i-K'2p5y 

:= °'<'-^'' (28) 

Also x'"^ = SVC? ^ = cos^ 6. Equation (27) determines the exact amount that 
the angular integration satisfies after a complete radial oscillation for a bound 
non-spherical polar orbit. 

Now using the latitude variable ^' the change in latitude after a complete 
radial oscillation leads to the following exact expression for the periastron ad- 
vance for a test particle in a non-spherical polar Kerr orbit, assuming a vanishing 
cosmological constant, 

^^GTR ^ A^- - 27r 

1 1 ^^^11. ..2\ a\l-E^)\ 



rrz J- J- " ^ / 



where the Abel-Jacobi's amplitude am(/u, k'^) is the function that inverts the 
elliptic integral 

= M (30) 



\/l - k'2 sin^ \1> 

in other words = amip, k'"^). For those astrophysical applications for which 

the modulus ^'^ = 0, A*^™ = ^^^J^^p (i, 1, 1, 

27r . If we assume an initial value for the latitude variable ^ different from 
zero then 

/ /7T 4 1 1 TT /I 1 

y VT^^Vc^VW^^ V2 2' 

. , ^ /I 1 1 3 . 2 , ,2 • 2 T \ a^(l-^^^ 



2' 2' 2' 2' 7 ' g / 

(31) 

and A*G™ = _ 27r = (* - *o) - 27r 

The orbital solution of non-polar non-spherical geodesies is discussed in Ap- 

pciidix C. 



^^Alternatively this can be seen from the fact that am{u, K?) = ArcSin(sn(M, k^)) and 
sn(u, 0) = sin(u). 

"/o*" ^i_^,f,i„2(^) = sin*o-Fi (|,^,|,f,sin2Vo,K'2sin2*o) see Appendix A.0.2 for 
the exact derivation of this angular integration. 
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2.1 Calculation of the orbital period in non-spherical po- 
lar Kerr geodesies 

In this subsection we shall perform an exact calculation of the period P for 
the non-spherical polar orbit of a test particle in Kerr gravitational field. The 
relevant differential equation is 

and we integrate from periapsis to apoapsis and back to periapsis. 

The physical quantity P, is an important rclativistic parameter, since orbits 
with small period and close to the event horizon exhibit large rclativistic effects 
which are easily measurable. The resulting expression has an elegant analytic 
form in terms of generalized hypergeometric functions of Appell and Gaufi's 
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hypergeometric function. We obtain 



ct = cP = 



2o GMbh 



r^(i/2) 



2'^' 2' 



2a;7r(3/2)r(l/2) 
^5 f(2^ 



,31 , 
Fi ( -,2,-,2,a;,«2 



7V r(5/2)r(i/2) 



r(3) 



0,2,i 3,w,k2^ 



+ £;a^2^ 



1 



+ 



'IEGMbh / w 



1-E^ ^{a~P){P-5) 



1.1. 

2>1> 2'!''^' 



a;7 



- ^-P^i U>l.o>2,a;,K' 



+ 



+ 



+ 



4:EG2Mbh 



AE2GMBii 



r(l)r(^) 
r(2) 

1 



u;3/2^ /3 1 22 



r(i)r(^) 
r(2) 



■^1 ( ^,1,^,1,«+'M' 



r(i) ^ 

r(i)rQ] 
r(2) 



a^E^^ 



1 _ 1 



2,1, 2' l''^-'^^ 



r(i) 



+ 



+ 



■ sin((^)Fi 



sin(<y9)Fi I i, i, i, ^,sin2v5,K2'sinV 



1113.2 2/ . 2 

sm (p, K sm 



2' 2' 2' 2 

• / X A 1 13.2 2/ . 2 M 1 
+ sin(v?)Fi (-,-,--,-, sin v'.K sin ) }> x — 



where k'^ is given by (28) and 

4 



V? = am I \fQ 



1 



1 TT „ /I 1 



V Vl - 2 V2'2 



1,K' 



2 \ ,.'2 



Also 



- 2r' 



rl — 



-a^ + 2r'_ 



r"_ — rj 



(34) 
(35) 



(33) 
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and H± := ^ {1 - E'^){a - f3){r'^ - I3)^JJ^. The moduli (variables) of the 
hypergeometric function of Appell are given by 



.2 



II = K 



2 ^ Q-/3(5-7 
a — 7 (5 — /3 
a — f3 r'j. — 7 



a — 7 rj_ — /? 

(36) 

The Lense-Thirring Period (LTP), is defined in terms of Lcnsc-Thirring's 
precession equation (12) and the orbital period P, equation (33), as follows 

LTP := -^TK (37) 

therefore its expression in closed analytic form is given entirely in terms of the 
generalized hypergeometric function of Appell Fi and Gaufi's hypergeometric 
function F. It constitutes the first exact expression for this important parameter 
for non-spherical polar Kerr orbits. 



2.2 Frame dragging of observed stellar orbits of S-stars 

We wish now to apply our exact formula, equation (12) for the calculation of 

the gravitomagnctic effect (Lcnsc-Thirring effect) for the orbits of S-stars in the 
central arcsecond of Milky Way assuming that the galactic centre region Sgr A* 
is a rotating Kerr black hole. 

We note that the careful study and observation of three-dimensional stellar 
orbits of S-stars at the central arcsecond of Milky Way provided the first strong 
evidence for the existence of a supermassive black hole at the galactic centre 
region Sgr A* [13, 14, 15, 16, 18]. Indeed, fitting of Kcplcrian orbits with 
experimental data revealed that the focus of the ellipse is essentially located 
at the radio position of Sgr A* As a matter of fact, Schodel et al [13] and 
Ghez et al [14] have shown that the focus of the orbit of the S-star S2 is located 
within a few milli-arcseconds of the nominal radio position of Sgr A* . The 
observed S-stars are B main-sequence stars as can be inferred from their H I 
Br7 absorption In particular S2 whose complete orbit has been observed 
appears to be a 15-20 Mq main sequence 08-BO dwarf star with an estimated 
age tg^2 < lOMyr [14] . A remark is in action: since ms <C Mbh i-e. the masses 

^*The J2000 equatorial coordinates of Sgr A* measured at 1996.25 were found to be at: 
R.A.=17''45™'°40^0409, decl.= -29°00'28".118 [21]. 
^^As is explained in [22], the observed stellar orbits can lead to a more precise determination 

of the distance from the Sun to the centre of the galaxy. It is well known that the distance 
_Ro from the Sun to the galactic centre not only is important for determining the structure of 
Milky Way but also holds an important role in establishing the extra galactic distance scale. 
It also affects estimates of the dark matter in the Local Group by affecting the Andromeda 
infall speed [23] . 

^^The wavelength of the Brackett-7 line is A(Br7 — line) = 2165. 5nm. 
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Q 


E 




LTP(yr) 


P (yr) 


SI 


71705.3106 


0.999993919 


0.0702036 ^'""l. 

revolution 


1.73 X 10^ 
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S2 


5693.30424 


0.999979485 


3.14297 T^ 

revolution 


6.25 X 10*5 


15.15 


S8 


9238.88965 


0.999992385 


1 51944 

revolution 


5.71 X 10^ 


67.0 


S12 


10641.9649 


0.999991213 


1.22889 T^ 

revolution 


5.70 X 10^ 


54.05 


S13 


36897.1620 


0.999988560 


0.190214 ''''=1. 

revolution 


2.48 X 10^ 


36.4 


S14 


5321.06355 


0.999988863 


3.47912 '^^.'=1 

revnliit.inn 


1.41 X 10^ 


37.9 



Table 1: Predictions for Lcnsc-Thirring effect Acj) , LTP and orbital period P 
for the observed orbits of S-stars around the galactic centre rotating supermas- 
sive black hole. The Kerr parameter was chosen to be acaiactic = 0.52^^^^. The 

values of Carter's constant Q are in units of — . 



star 


Q 


E 




LTP(yr) 


P(yr) 


SI 


71705.3106 


0.999993919 


0.134183 T!- 

revolution 


9.07 X 10» 


93.9 


S2 


5693.30424 


0.999979485 


6.0073 

revolution 


3.27 X 10*^ 


15.15 


S8 


9238.88965 


0.999992385 


2.90417 

revolution 


2.99 X 10^ 


67.0 


S12 


10641.9649 


0.999991213 


2-34883 


2.98 X 10^ 


54.05 


S13 


36897.1620 


0.999988560 


0-363565 ..^r^^ 


1.3 X 10^ 


36.4 


S14 


5321.06355 


0.999988863 


6.64981 ^--iL. 


7.39 X 10^ 


37.9 



Table 2: Predictions for Lense-Thirring effect Acj) , LTP and orbital period 
P for the observed orbits of S-stars around the rotating galactic centre black 
hole. The Kerr parameter -was chosen to be acaiactic = 0.9939^^. The values 
of Carter's constant Q are in units of ^-^^sn . 



of S-stars are very small compared to the central black hole mass the effect 
of their gravitational fields on the black hole's gravitational field is negligible. 
Therefore the approximation that they move on a timelike geodesic of the four 
dimensional world is entirely justifiable and consistent . They can be treated as 
test particles. 

We now proceed to calculate, using our exact analytic solutions, the Lense- 
Thirring effect and the corresponding Lense-Thirring period for the observed 
orbits of S-stars: S1,S2,S8,S12, S13,S14 for various values of the Kerr parameter. 
These orbits have been measured by Eisenhauer et al [18] at near-infrared, with 
the data taken with the new adaptive optics-assisted integral-field spectrometer 
SINFONI on the European Southern Observatory (ESO) Very-Large- Telescope 
(VLT). Orbital data for the semi-major axis a and the eccentricity e of the corre- 
sponding Keplerian orbits are listed in [18]. Our motivation for this application 
was to investigate if the dragging of the orbital plane by the Lense-Thirring 
effect can help in understanding the observed distribution of the orbital planes 
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star 


Q 


E 




LTP(yr) 


P(yr) 


SI 
S2 
S8 
S12 
S13 
S14 


71705.3106 
5693.30424 
9238.88965 
10641.9649 
36897.1620 
5321.06355 


0.999993919 
0.999979485 
0.999992385 
0.999991213 
0.999988560 
0.999988863 


0-134372 
6.01576 ^f=^ 

revolution 

2.90826 T^ 

revolution 

2.35213 T^ 

revolution 

0.364077 

revolution 

6.65917 Tt- 

revolution 


9.05 X 10** 
3.26 X 10^ 
2.98 X 10^ 
2.978 X 10^ 
1.295 X 10^ 
7.37 X 10^ 


93.9 
15.15 
67.0 
54.05 
36.4 
37.9 


Table 3: Predictions for Lense-Thirring effect A(j)'~^'^^, LTP and orbital period 
P for the observed orbits of S-stars around the rotating galactic centre black 
hole. The Kerr parameter was chosen to be acaiactic = 0.9953^^^1. xhe values 

of Carter's constant Q are in units of — • 


star 


Q 


E 




LTP(yr) 


P(yr) 


SI 

S2 

SB 

S12 

S13 

S14 


71705.3106 

5693.30424 
9238.88965 
10641.9649 
36897.1620 
5321.06355 


0.999993919 

0.999979485 
0.999992385 
0.999991213 
0.999988560 
0.999988863 


0.134978 ^[^^ 

revolution 

6-042901 j,,^ 
2.92138 T!- 

revolution 

2.36275 "'^^ 

revolution 

0.36572 T^ 

revolution 

6.68921 '^''=1. 


9.01 X 10** 
3.25 X 10^ 
2.97 X 10^ 
2.96 X 10^ 
1.29 X 10^ 
7.34 X 10*^ 


93.9 
15.15 
67.0 
54.05 
36.4 
37.9 



Table 4: Predictions for Lense-Thirring effect , LTP and orbital period P 

for the observed orbits of S-stars around the rotating galactic centre black hole. 
The Kerr parameter was chosen to be acaiactic = 0.99979 ^^°" . The values of 

Carter's constant Q are in units of — 33™ • 



star 


Q 


E 


A0^'^'«' 


LTP(yr) 


P(yr) 


SI 

S2 

S8 

S12 

S13 

S14 


65431.43821 
5273.53220 
6524.03518 
9730.48589 
26279.90639 
4204.76359 


0.999993545 
0.999979145 
0.999991945 
0.999990853 
0.99998444 
0.999987653 


0-15394,,,^ 
6.73964 

revolution 

4.89646 T^ 

revolution 

2.68671 Tf- 

revolution 

0-604888^;^^ 
9.47151 Tt- 


7.22 X 10* 
2.84 X 10*^ 
1.63 X 10^ 
2.45 X 10^ 
4.91 X 10^ 
4.43 X 10'^ 


85.8 
14.8 
61.6 
50.9 
22.9 
32.4 



Table 5: Predictions for Lense-Thirring effect A(^ , LTP and orbital period 
P for the observed orbits of S-stars around the rotating galactic centre black 
hole. The Kerr parameter was chosen to be acaiactic — 

0.9939^^. The values 

of Carter's constant Q are in units of ^-^^sn . 
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for the S-stars. The parameters of the exact theory of non-spherical polar orbits 
developed in the previous sections are the invariant parameters Q,E, the spin 
and the mass of the galactic black hole, a, Mbh respectively. For the black hole 
mass wc adopt a value: Mbh = 4.06 x IO^Mq. The choice of values for the 
invariant parameters is restricted by requiring that the predictions of the theory 
for the periastron, apoastron distances and the period P are in agreement with 
the orbital data in [18]. We present our results for four different values of the 
Kerr parameter in Tables 1-4 respectively . Indeed as we can see the stars 
S2 and S14 have LTP's of the right size as compared with estimations of their 
age from spectroscopic analysis and the average stellar ages in the two outer 
star rings/discs [19]. According to [19] it is estimated that the stellar discs in 
the galactic centre have formed about 6 ± 2 Myr ago. Thus S2 and S14 be- 
cause their LTPs are comparable to the average ages of the stellar contents of 
the discs could have originated from the stellar discs before the Lense-Thirring 
dragging brought them where they are currently observed. The LTP's decrease 
with higher values of the spin of the black hole. Our results in Table 1 agree 
with the results in [18], where the Lense-Thirring periods were calculated for 
a = 0.52 using a weak field approximation formula [20]. This is a nice check of 
the weak field limit of the exact theory. 

We repeated the analysis for a = 0.9939 but for different values of the 
invariant parameters Q,E. This choice is still consistent with experimental 
data for the semi-major axis a and eccentricity e quoted in [18]. The results 
are presented in Table 5. Now we see that also the LTP for S8 is of the right 
magnitude. Thus we can conclude that the Lense-Thirring effect can help in 
understanding the observed distribution of the orbital planes for some of the S- 
stars but perhaps for not all of them On the other hand when more precise 
astrometric data will become available for all of the S-stars one should try to 
fit the experimental raw data using the relativistic geodesic equations of motion 
and the exact results obtained in this work, assuming that the galactic centre 
is a Kerr black hole, rather than try to fit the data with Keplerian orbits. 

Our choice of value for the black hole mass is consistent with the best fit 
central mass (for a distance Rq to the galactic centre of 8Kpc) in [18]. If we try 
to fit the data for the periastron and apoastron distances in [18] with a smaller 
black hole mass we find that the orbital period P increases. In particular, 
adopting the value Mbh = 3.59 x IO^Mq reproducing the measured values for 
rp,rA [18] for the star S2 results in P> 15.6yr in conflict with experiment [18]. 

^''^In tables 3,4 we have calculated the Lense-Thirring precessions for the values of the spin of 
the black hole for which the Aschenbach effect (hump of the LNRF-velocity profile) becomes 
to be relevant a=0.9953 for Keplerian discs [27] and a=0. 99979 for marginalUy stable discs 
respectively [28]. 

^®The periapsis advance for the S-staxs with the initial conditions as those in tables 1, 2 
have been calculated, using the exact closed formula eq.(29), in appendix E, table 19. 
^^See [19] for a recent summary of attempts to explain the youth of S-stars. 
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parameters 


Periastron 


Apoastron 


predicted dragging 


aCalactic = 0.9939, 


5.51905 


14.5373 


A</>GTR ^ 108552 = 223905^5^^^ 


E = 0.956853, Q = 13 






=62.19°per revolution 


ClGalactic = 0-52, 


4.80864 


14.5620 


^ ,GTR _ i Q^^g ^ 213464 '^'Z 

" revolution 


E = 0.956853, g = 13 






=59.29'^per revolution 



Table 6: Predictions for frame dragging A<j> from galactic black hole for 
a timelike non-spherical polar stellar orbit for two different values of the Kerr 
parameter acaiactic- Carter's constant is fixed at Q = 13 — jjSh -^^iiiie E = 
0.956853. The values of the radii and Kerr parameter are in units of GMbh/c^- 

3 Zooming in on the galactic centre black hole 

3.0.1 Frame dragging and periapsis advance of non-spherical polar 
orbits in the central milliarcsecond of Sgr A* 

In this subsection we shall apply our exact formulae for the calculation of the 
relativistic periastron advance and frame dragging for non-spherical polar orbits 
in the central milliarcsecond of the galactic centre of Milky Way. These cusp 
orbits are the target of the GRAVITY experiment and will provide an interesting 
enviroment for probing stellar dynamics in the strong-field regime of General 
Relativity ^o. 

Assuming that the centre of Milky Way is a rotating black hole and the 
spacetime near the region Sgr A* is described by the Kerr geometry, we cal- 
culated using (12) the precise frame dragging of a timclikc stellar orbit with a 
non-spherical polar geometry. The results are displayed in tables 6 and 7. The 
periods for the orbits in table 7 were calculated to be: 6.443 x 10^ s, 6.359 x 10^ 
s for a = 0.52 '^^bh ^ 9 9939 gMbh respectively, assuming a galactic black hole 
mass Mbh = 4.060 x IO^Mq. " 

Let us give subsequently an example of a non-spherical stellar polar orbit of 
low eccentricity. The results are exhibited in table 8. 

Using the exact formula, equation (29), the periastron advance for the or- 
bits with initial conditions displayed in table 8 were calculated. Our results 
arc summarized in table 9. Thus, the closer we get to the galactic centre 
the relativistic effect of periapsis advance becomes very substantial. We ob- 
serve that the effect increases with decreasing Kerr parameter. The semi-major 
axis a = ^^'^'^^ ~ HQ^I^pji = 55i?,, corresponds to an angle of observation 
~ 0.54TOas = 540/Ltas (lmas= 10"'^ arcs)-^^. The calculated orbital period P is 
1.72 days. Measurement of such orbits and the corresponding relativistic effects 

^"GRAVITY is an adaptive optics (AO) assisted, near-infrared VLTI instrument for preci- 
sion narrow-angle astrometry and interferometric phase referenced imaging of faint objects. 
With an accuracy of 10 micro-arcseconds (10/^as), GRAVITY will be able to study motions 
to within a few times the event horizon size of the Galactic centre black hole. 

^^We assume a black hole mass Mbh = 4.06 x IO^Mq and a galactic centre distance of 
S.OKpc. 
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parameters 


Periastron 


Apoastron 


predicted dragging 


aCalactic = 0.9939, 


7.90470 


12.5692 


A0«"' = 0.6477498 = 133608^^^^ 


E = 0.956853, Q = 14 






=37.11°per revolution 


aCalactic = 0.52, 


7.66620 


12.6339 


A0°™ = 0.366361 - 75567.5 

^ revolution 


E = 0.956853, Q = 14 






=20.991°per revolution 



Table 7: Predictions for frame dragging from galactic black hole for timelike 
non-spherical polar stellar orbit for the same values for the Kerr parameter and 

E as in table 6 but for a value of Carter's constant Q = 14 — • The values 
of the radii and Kerr parameter are in units of GMbh/c^- 





parameters 


Periastron 


Apoastron 




predicted dragging 




aCaiactic = 0.99979, 


107.703 


112.222 


= 0.0112114 - 2312.52 

revolution 


E 


= 0.9954853, g= 113 








=0.642368°per revolution 




dGalactic = 0.99616, 


107.703 


112.222 




= 0.0111707= 2304.13 T^ 

revolution 


E 


= 0.9954853,0== 113 








:0.640036°per revolution 




aCalactic = 0.9953, 


107.703 


112.222 


A<^«TR 


= 0.0111611 = 2302.14 

revolution 


E 


= 0.9954853,(5= 113 








=0.639484°per revolution 




aCalactic = 0.9939, 


107.703 


112.222 


A</,G™ 


= 0.0111454= 2298.9 T!- 

revolution 


E 


= 0.9954853,0= 113 








=0.638584°per revolution 




OGalactic = 0.52, 


107.697 


112.227 


A(/.«™ = 


= 0.00583177= 1202.89^5^;^ 


E 


= 0.9954853,0= 113 








=0.334136°per revolution 



Table 8: Predictions for frame dragging from a galactic black hole for timclikc 
non-spherical polar stellar low- eccentricity orbits. The values of the radii and 
Kerr parameter are in units of GMbh/c^, while those of Carter's constant Q in 
units of (GMbh/c2)2. 
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pareimeters 

OGalactic = 0.99979, 

= 0.9954853, = 113 
aCalactic = 0.99616, 
= 0.9954853, g= 113 

OGalactic = 0.9953, 

= 0.9954853, Q = 113 
aCalactic = 0.9939, 
i; = 0.9954853,(5= 113 

flGalactic = 0.52, 

= 0.9954853, Q = 113 



rp 
107.703 

107.703 

107.703 

107.703 

107.697 



112.222 
112.222 
112.222 
112.222 
112.227 



predicted periastron advance 
^^GTR ^ 0.178388 = 36795.2 



revolution 

=10.2209°per revolution 
A^GTR ^ 0.178391 = 36795.9 T!- 



revolution 

=10.2211°per revolution 
^^GTR ^ 0.178392 = 36796.0 T!- 

revolution 

= 10.2211°per revolution 
A*c:tr ^ 0.178393 = 36796.3— a^f^ 

revolution 

=10.2212°per revolution 
A*GTR ^ 0.178724 = 36864.4 ""'"l 

revolution 

=10.2401°per revolution 



Table 9: Predictions for periastron advance from a galactic black hole for time- 
like non-spherical, polar, stellar low eccentricity orbit. The values of the radii 
and Kerr parameter are in units of GMbh/c^, while those of Carter's constant 
Q in units of (GMbh/c^)^. 

close to the outer event horizon r+ of the candidate galactic black hole can 
provide an important test of the theory of general relativity at the strong field 
regime. 

A very interesting example of a stellar polar non-spherical orbit is presented 
in Table 10. We present the orbital parameters for two choices of the spin 
of the black hole. The rest of the initial conditions are g = 380 "^ 4"" ,^ = 
0.99954853. The angular resolution of the pcriapsis rp is ~ Imas. We see from 
the results in the table that the relativistic phenomenon of pcriapsis precession 
is quite substantial and it amounts to with the orbital period of only 53.6 
days. For the convenience of the reader we also list the pcriapsis precession and 
Lcnso-Thirring effect after one revolution as well as the Lense-Thirring period. 
As we mentioned, VLBI observations aim to detect such cusp stellar orbits which 
are closer to the candidate black hole than the observed S-stars. The pcriapsis 
precession in the example we discuss should be measured with a good precision 
after one year of observations. We repeated the analysis for higher values of the 
spin of the black hole and the results are dislayed in table 11. 

One might also wonder about the tidal effects on the orbits of the stars we 
just discussed. We can estimate these effects as follows: tidal disruption of the 
star with mass m* and radius i?* is expected to occur only if the star approaches 
the black hole to within its Roche radius 

rn=(^Y'R. (38) 

For a star like our Sun with i?* = Rq and = Mq and for a galactic black 
hole mass Mbh = 4.06 x IO^Mq the Roche radius is calculated to be = 
18.5 *^^°" ■ Thus the ratio of the periastron radius to the Roche radius for the 
orbits presented in tables 9,11 is 5.8 and 11.2 respectively. Consequently, any 



19 



Orbital parameters 


Case 1: a=0.52 


Case 2: a=0.99616 


rp 


207.628 


207.628 




2005.837 


2005.837 


Periapsis advance(arcs/rev.) 


10460.03 


10454.7'' 


Frame dragging (arcs/rev) 


186.626 


357.513 


Periapsis advance (^) Assumes Mbh = 4.06 x lO^M© 
L.T. effect (^) 
Orbital Period P (days) 


19.8 


19.8 


0.353 


0.675 


53.6 


53.6 


L.T. Period (yr) 


1.02 X 10^ 


533 



Table 10: Calculation of orbital data and the corresponding relativistic effects 
for a cusp polar non-spherical stellar orbit, around the galactic black hole, for 
two different values of the Kerr parameter with Q = 380^^^, E = 0.99954853. 



Orbital parameters 

rp 

Periapsis advance(arcs/rev.) 
Frame dragging (arcs/rev) 
Periapsis advance 



L.T. effect (^) 
Orbital Period P (days) 
L.T. Period (yr) 



^) " Assumes Mbh 



4.06 X 10'=' Mp 







a=0.9953 

207.628 
2005.837 



10454.8 
357.204 

19.8 
0.675 

53.6 

533 



a=0.99979 

207.628 
2005.837 
10454.7'' 
358.815 
19.8 
0.679 
53.6 
531 



Table 11: Calculation of orbital data and the corresponding relativistic effects 
for a cusp polar non-spherical stellar orbit, around the galactic black hole, for 
two different values of the Kerr parameter with Q = 380 ' i'^" , E = 0.99954853. 
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tidal effects are small for a solar-type star in orbit around the galactic centre 
black hole according to the data in tables 9, 11 and we do not expect significant 
mass loss. If a star similar to S2 with m, = 15Mq, i?* = 7Rq is in a polar 
orbit in the central mas of Sgr A* according to the data in tables 9, 11 the 
ratio ^ is calculated to be 2.05 and 3.95 respectively. For a neutron star with 
mass ~ IAMq and radius ~ lOKm in a non-spherical polar trajectory with the 
orbital characteristics of tables 9-11 the tidal effects arc negligible. 

On the other hand in order to investigate more precisely the tidal interaction 
of a star in a high eccentricity orbit around the supermassive galactic centre 
black hole a full relativistic treatment is required. Such an analysis which will 
involve the calculation of the tidal tensor for non-spherical polar, equatorial 
or arbitrary orbital inclination test particle trajectories with or without the 
cosmological constant is however beyond the scope of the ciirrent work. It is 
a promising avenue for future research since the tidal interaction of a star by 
supermassive black hole may lead to observable effects [36]. 



4 Periastron precession of non-circular equato- 
rial orbits around a central rotating mass in 
the presence of a cosmological constant 

We now proceed to derive exact expressions for the periastron advance and 

orbital period for a timelike non-circular orbit in the Kerr field in the presence 
of a cosmological constant. In the former case the relevant differential equation 

is 22 

dr ArVR ^ J ^ ^ 

In (39) R is defined in (5) with Q = 0. 

For the derivation of an analytical exact expression of periastron advance 
for the equatorial orbit of a test particle in Kerr-de Sitter spacetime we need to 
calculate the following Abelian integrals 



r -^i^^-J^dr + r ^=5= \{r^ + a^)E - aL] dr (40) 



Let us now calculate the above Abelian integrals. The first integral is in the 
general case (when all roots are distinct) a hyperelliptic integral. By applying 
the transformation 



and organizing the roots of the polynomial R as follows 



(41) 



> > Q!p > ai (42) 



^^Setting Q = 0,6 = n/2 in the relevant equations in equations 5. 
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where = a^^i^Up = a^+i = rp,a^ = ta^oli = a^+j+i,? = 1,2,3, our genus 
2 hyperelliptic integral is reduced to the integral representation of Lauricella's 
fourth, generalized hypergeometric function of three variables Fjy [37] More 
precisely we obtain 

+ 



-Fd 



1111 



2' 2' 2' 2 

r(i)r(i) 



, 1, K , A , fj," 



r(i_ 
r( 



r(2) 



where 

H := (a^ - a^+i)(a^ - a^+2)(a^ - a^+3)(a^ - a^+i) (44) 
and a; is now defined as 

w := "'^ ~ °'^+' (45) 

In addition we have defined the dimensionless quantity A from the relation 
A :=A( '^^^" )^. The variables of the hypergeometric function Fd are given by 
the expressions 

Q!/i-i - a^+i - ol^ - a^+2 "r)^- ^ a - 7 

,2 Q^M-i ^ a - /? ''a 

A — UJ- 



- a/i+3 - /3 a 

2 Q - /3 4 - 
M - 7 Ti 3-:^ — -:a 

A 



(46) 



Now for the second integral in (40) we use partial fraction expansion 



h2 [(r^ + a2)F - ail A^ 

^ ^ - :r + 3 + + (47) 



r — r J r — r^^^ r — r+ r — r. 



^^See Appendix B for its definition, the differential equations it obeys and its integral 
representation. 
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For the coefficients 
A' 



A^ = 



1, • • • 4 we obtain the expressions 

-3a2s2 [aE-L + Eaf] 

(tti - - 7i)(«i - Si)A 

Sa^S^ [aE-L + Edf] 
(ai-(5i)(Ji-/3i)(,5i-7i)A 

Sa^S^ [aE-L + Efjj] 
(-ai+/3i)(-/3i+7i)(/3i-(5i)A 

-Sa^S^ [aE-L + E^j] 
(ai-7i)(/3i-7i)(7i-(5i)A 



(48) 



where we have defined 



Q!i,r-+ = /3i,r_ = 7i,rA = (5i 



(49) 



for the roots of the polynomial A,-. Proceeding in a similar way as before 

we transform the Abclian integrals into the familiar integral representation of 
Lauricella's function Fd ■ For instance for the integral J^^ ^^^^ ^ we derive 
an exact expression in terms of Lauricella's function Fn of four variables: 



A4 



dr 



A, 



(r. 



+ 



- 2uFd 



2' 2' ' 2' 2 ''^ 

1 1 i i 9 K-2 x2 2 2 

2' 2 2' 2 



r(i)rQ) 
r(3) 



r(i)r(i) 



r(2) 



(50) 



In equation (50) the four moduli are given in terms of the roots of the polynomial 
R and the radii of the black hole horizons by the expressions 



an 



a,. 



-1 — Oiij,+2 



a,, 



'A 

a - 



-/3 7- 

(3 r_ 



a 



a 



fj- ' A 
a-(3 rl 



a^_i - a^+i an+6 - a^ 



f3 r 



(51) 
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while H- is defined to be: 



= ^(a - /3)(a - 7)(a - r-)Ja{a - ) (52) 
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In total we obtain: 

J-E^iaE-L) 



A(l)' 



GTR 

E 



+ 



+ 



+ 



+ 



+ 



+ 



F /'I i i i 1 " 
\ 2, 2' 2' 2' ' 



/3 



A 



/3 rl 



3 111 



2'2'2'2'^'ri 



r\ — P a — — /? a 

a — /3 rj^ — 7 a — /? r 



— /3 a 

/3 4 - 



A 



/3 a — 7 



/3 a 



/3 a 



S-'^A 



2uFd 



5.111 ^ 7a-^J 
2'^'2'2-2'^'^7l^ 

3 1 i i i 2 ^ ^A-"A 

2 2 2 2 rT — a 



7 — a 



— ,UJ— 

a a 



'1 111 r — r 



7 — a 



,L0- 



a 



r(i/2)^ 
r(i) 

r(|)r(^- 
r(2) 



r(i)r(|) 
r(3) 

r(i)rQ) 
r(2) 

r(i) 



t-f^A 

2 /5 1 1 1 7- 
'"^^U'2'2'2'^'^'^T 



a 



A „ TA-^n r(|)r(i) 

r(l)r(i) 



— ,CJ — — 

Q a r1 



1111 



7- 



2' 2' 2' 2'^'^''^7 



— ,UJ — 

a a 



„2 



,UJ- 



r(2) 

r(^)r(^) 
r(i) 



3^+ 



2'2' '2'2"^''^7-a 



, a; — , w— 9- 

r+ — a a 



r(3) 



/3 1 11 7- 
2'^^-l2'2'^'2'2'2'"T 



-,u;- 



r+ — a 
^1 



-, a; — ,u)- 



a 



a 



p fi 1 1 1 1 1 r+ 

\22 22 7 — a r+ — a 



2 1 
A 'A ''a ^ A 



n 

a 



n - a 



r(2) 

rQ)r(i) 
r(i) 



/5 1 1 1 

Fd I 7:>7:'1> 2' 2' ' 



7- 



r 



- 2a;FD 



^2' 2 

3 111 7 
2- 2'!' 2' 2'^'"^ 7^ 



a 



-,UJ- 



1 2 
A , .''a , ,^A 

q; a 



A 



r(l) 



r(3) 
r(i)r(i: 



,'11 11 7-ri 
+ ^-(2'2'S'2'1'^^ 




,0; 



where the Hj^, Hj^ are given by the expressions 



Hi = (a-r+)^(a-/?)(a-7Ha-ri) 



(54) 



The periastron advance is then given by the expression 



A<^g™ - 2w (55) 

The phenomenological apphcations of equations (53) and (55) will be performed 
systematically in a future publication 

4.1 Exact, closed form formula, for the orbital period of 
equatorial non-circular Kerr geodesies 

In this subsection we assume vanishing cosmological constant. For equatorial 
orbits Q vanishes and 6 = 7r/2. The relevant differential equation for the calcu- 
lation of the orbital period Pe for equatorial non-circular Kerr geodesies is 

dt _ jr'^ + a'^)[E{r^ + a"^) - aL] a{L - aE) 
Then we want to calculate exactly the following definite integral 

In (56), (57) the qiiartic polynomial R is defined in (5) with A= 0, Q = 0. The 
idea is to use the appropriate transformation and bring the above radial integral 
into the integral representation of generalized hypergeometric function of Appell 



*The calculation for the corresponding orbital period is outlined in appendix D. 
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Fi. We then obtain ^5 



ct = cPe 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



2cJ7r(3/2)r(l /2) 
^5 W) 



2' '2 



T^c.^ r(5/2)r(i/2) 



2E2GMbh 



f3 



1 



l-£;V("-/9)(/3-'5) 

r(i)r(i) 



r(i) 



r(3) 



1 1 



Fi 



4S2GMBH 
4E2GMBU 



r(2) 
1 



3/2^_ 



^(a-/?)(/3_5) 1^2' 2' ' y r(i) 

H+ 
1 ^ 1 



r(i)rG) 
r(2) 



-^2 /l^ 



+ 



i7_ 



1 T 1 1 2 2 



r(i) 

r(l)r(|) 
r(2) 

r(i) 



+ aL2 



GMbh 



^3/2^', /3 1 2 2 

^Fif-,l,-,2,4,/.2 



r(i)r(^) 
r(2) 



-2 /l^ 



+ 



-Fi 



1 1 1 2 2 



r(i) 

r(f)r(£ 
r(2) 

£l(i) 
r(i) 



where 

:= ^- ^= (59) 

and A± are given by equation (35). Also rj_ := 1 ± Vl — a2, a arc the dimen- 
sionless horizon radii and spin of the black hole respectively. Equation (58) is 



^^See appendix A for some details of the calculation. 
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the first exact expression in closed analytic form of the period of a test particle 
in a non-circular equatorial orbit around a Kerr black hole in terms of Appell's 
generalized hypergeometric function Fi and Gaufi's ordinary hypergeometric 
function. It constitutes a generalization of the Kcplerian period for a circular 
equatorial orbit of radius r in the Kerr field [35] (see also Eq.(105) [11] for A = 0) 

a /iMIm ± ,3/2 
t = 2Tr—^ ' (60) 

Therefore equation (58) is of paramount importance for calculating the orbital 
and periapsis precession periods for non-zero eccentricity equatorial orbits of 
test particles in the vicinity of the horizon r+ of the Kerr black hole in the 
strong-field regime of general relativity. In addition, it may be of importance 
for the phenomenology of quasi periodic oscillation (QPO's) frequencies [39, 27] 
and accretion physics. 



4.1.1 Examples of equatorial non-circular orbits with short orbital 
period 

In this subsection we calculate the periapsis advance for equatorial non-circular 
orbits with short orbital periods in the asymptotically flat Kerr black hole. In 
this case the contribution of the spin of the black hole to the relativistic effect of 
periapsis advance is maximized as compares to the polar orbits. Consequently, 
a precise measurement of the periapsis advance can lead to an independent 
measurement of the spin of the black hole. The exact expression for the periapsis 
advance for the orbit of a test particle in a non-circular equatorial motion around 
a Kerr black hole has been derived in [12] and it involves the hypergeometric 
function Fi of Appell By applying this formula for the choices of initial 
conditions L = ^/380^^^,E = 0.99954853 and for two values of the Kerr 
parameter we obtain the results displayed in Table 12. For the calculation of 
the equatorial orbital period Pe, eq.(58), for this high-eccentricity equatorial 
orbit we assumed a galactic black hole mass Mbh = 4.06 x IO^Mq. 

Indeed we observe now the important contribution of the Kerr parameter on 
the relativistic effect of periapsis advance. 

A more dramatic example of an equatorial non-circular orbit closer to the the 
event horizon r+ of the galactic rotating black hole is presented in table 13. From 
the table we see orbital periods in the range 105min-108min depending on the 
value of Kerr parameter and quite substantial relativistic periapsis advance. The 
calculated orbital frequencies are: 0.155mHz, 0.158mHz for a = 0.52,0.99616 
respectively. For comparison we repeated the calculation for a microquasar with 
black hole mass Mbh = lOM© and for a = 0.99616^^, E = 0.9586853, i = 
^^/j^ GMBH _ 'j'jjg calculated orbital frequency has the value 64.2Hz. 

^®For the convenience of the reader we exhibit this exact formula in appendix E. 
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Orbital parameters 


Case 1: a=0.52 


Case 2: a=0.99616 


rp 


207.748 


207.855^^ 




2005.825 ^^F^ 


2005.813 


Periapsis advanceiAt/)^'^^ — 27r(arcs/rev.) 


IOO82.4'' 


9742.34^" 


Periapsis advance (^) Assumes Mbh = 4.06 x lO^M© 
Orbital Period Pe (days) 


19.05 


18.4 


53.6 


53.6 



Tabic 12: Calciilation of orbital data and the corresponding relativistic ef- 
fects for a cusp equatorial non-circular stellar orbit around the galactic black 
hole, for two different values of the Kerr parameter with L = \/380 '^^'^" , E = 
0.99954853. The exact expression A(?if ™ is given by Equation (111) in [12]. 



Orbital parameters 

rp 

TA 

Periapsis advance:A(/)E'^^ — 27r(arcs/rev.) 

Orbital Period Pe (s) Assumes Mbh = 4.060 x 10*^ Mq 

Orbital frequency z^E(mHz) 



Case 1: a=0.52 
7.49878 
15.1106 

497734"" 
6.462 X 10^ 
0.155 



Case 2: a; 
9.20674 



=0.99616 

GMbh 



14.1781 

309217^ 
6.319 X 10^ 
0.158 



Table 13: Calculation of orbital data and the corresponding relativistic effects 

for a cusp equatorial non-circular stellar orbit around the galactic black hole for 
two different values of the Kerr parameter with L = ^/\2^^^,E = 0.9586853. 
The exact expression At^^^^ is given by Equation (lll)in [12]. 



29 



5 Conclusions 



In this work we have investigated the motion of a test particle in Kerr spacetime 
with and without the cosmological constant. 

More specifically, wc have studied and derived for the first time exact so- 
lutions of the relativistic geodesic differential equations of motion for two im- 
portant classes of possible orbits: a) non-spherical polar Kerr orbits and b) 
non-circular equatorial Kerr orbits in the presence of the cosmological constant. 

In the former case the exact orbital solution is given in terms of the square 
of Jacobi's sinus amplitudinous elliptic function. Furthermore we have derived 
exact novel expressions for the Lense-Thirring effect, periapsis advance and 
orbital period. The resulting solution in closed analytic form for the Lense- 
Thirring gravitomagnetic effect has been expressed in terms of Appell's general- 
ized hypergeometric function Fi. The exact formula for the relativistic periap- 
sis advance has been expressed elegantly in terms of Abel's- Jacobi's amplitude 
function am{u, k^) [33] with the argument u given by Gaui3's hypergeometric 
function F and the initial conditions for the invariant parameters Q, E and the 
Kerr parameter. The closed form formulae for the orbit and the periapsis pre- 
cession have been achieved by using the idea of inversion of elliptic integrals 
[32] and properties of the modular functions. The derived solution for the polar 
orbital period P has been expressed in terms of Appell's generalized hypergeo- 
metric function Fx of two variables and GauB's hypergeometric function F. This 
allowed us to define the LTP (frequency) of polar non-spherical orbits around 
a Kerr black hole. We also derived the exact solution for the orbital period Pe 
for the equatorial non-circular orbit of a test particle in Kerr spacetime. The 
resulting formula again rests in terms of hypergeometric functions of Appell and 
Gaufi. 

Assuming that the galactic centre is a rotating Kerr black hole as is indicated 
by recent observations of the galactic centre, we applied our exact solutions for 
the precise frame dragging and periapsis advance of stellar non-spherical polar 
Kerr orbits near the outer event horizon of the supermassive black hole. More 
specifically, we investigated orbits in the central milliarcsecond of the galactic 
centre of Milky Way. Such orbits have short orbital periods and therefore if 
detected they can provide a unique laboratory for measuring with high precision 
the phenomenon of periapsis advance and frame dragging at the strong field 
realm of general relativity. Such angular resolutions will be achieved by the 
GRAVITY experiment using VLBI techniques. We provided examples with 
orbital periods in the range lOOmin-54 days. 

We have also calculated the frame dragging and periapsis advance effects 
for the more distant (with respect to the outer horizon of the Kerr galactic 
black hole) S-stars in the central arcsecond of Sgr A* for various values of the 
Kerr parameter and the other invariant parameters. Since the orbital periods 
in this case are in the range 14.8 yr-94 yr, more precise astrometric data from 
GRAVITY will allow the probing of relativistic effects on a longer (~ 10 year) 
time scale. 

In a similar fashion the exact solution of non-spherical non-polar orbits not 
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confined to the equatorial plane has been derived and was expressed in terms 
of Jacobi's sinus amplitudinous elliptic function. 

We have also calculated an exact expression for the periastron advance of 
a test particle in a non-circular equatorial orbit around a rotating mass in the 
presence of the cosmological constant (Kerr-de Sitter gravitational field). The 
resulting expression was given in terms of Lauricella's generalized hypergeo- 
metric functions F]j of three and four variables. This solution can be used for 
investigating in a precise way the combined effect of the cosmological constant 
and the rotation of the central mass on the orbits of test particles. Our exact 
solutions of the equations of general relativity involve important functions of 
mathematical analysis and the exploration of their topological aspects consti- 
tutes a very interesting topic for future research [42]. We find very pleasing 
the productive and exciting interplay of observational astronomy of the galactic 
centre with the exact theoretical strong-field regime results presented in this 
work. 
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A Precise calculation of the orbital period of 
non-circular equatorial geodesies 

Equation (57) can be written 

cPe = 2 E(r^+a^)^ — -^dr-2 — !^ — —=dr+2 

Jrp As/R Jrp As/R Jrp s/R Jrp 

(61) 



31 



Let us outline the calculation of the above definite integral. The first integral 
in equation (61) can be written 



/ ^ — 1= — dr 

Jrp s/R 



r2+a2-2GMBH:& 



rp y 11 Jrp y ii \ / Jrp 

(62) 

where 

r£(l!±f!)*+f (63, 



VR Jrp y/R c2 



Now the term 4£J ^ g^bh ^^ j^^ ^ using a dimensionless variable r'i.e. r = 
j.iG^H g^j^jj ^jjg transformation 

Lor' - a; r - 7 

with LV = ) is transformed into the integral representation of Gaufi's hyper- 
geometric function 



= 4E 



\ J Jrp VR 

GMbh dz 



c2 ^^i-E^){l3-5){a-P) Jo ^z{l-z){l-^iH) 

(65) 



c2 ^{l-E-^){(3-5){a-l3) V2' 2' ' 



Subsequently let us calculate 

drEf_ drEa^ drE ^GMbut 

Jrp s/R Jrp s/R Jrp s/R 



la 



The second term in equation (66) cancels the third term in equation (61). Now 
the first term in (66) using the transformation (64) can be transformed into the 
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integral representation of the generalized hypergeometric function Fi of Appell 
'■'■^ drEr"^ 

irp 



f 



R 

EGMbhV^ 



dz{'yivz — 



c^yJ{l-E-^){(i-5){a-p) Jo il-ujz)^y/z{l-z){l-fi^z) 



EGMbu\/(^ 



2 2 

7 w 



13' 



r(i)r(i) 



2w7 



1 



- ( o'2,o'2,w,A*^ 



r(i) 

r(i)rQ) 
r(2) 



(67) 



In a similar fashion the integral ^7= ^^'H""'' is calculated to be 



2£; 



,11 



£!(i) 
r(i) 



r(l)r(i) 



r(2) 



^F, (1 1,1, 2, a;,/. 
V(l--B2)(/J-^)(a-/3) 



where wc used the transformation (G4). 
Finally, the last term in equation (62) 



r(i) 



(68) 
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+ 



+ 



4E 



'GM] 



BH 



r-A n2 



-2GMbh: 



4EGMbh 



H+ 
H- 



Tp 

uf'I^A, 



-dr 



3 1 2 2 



+ 



-Fi 



£111) 
r(i) 

r(i)r(i) 
r(2) 

£!(i) 

r(i) 



r(i)r(i) 



r(2) 



(69) 



where H± := ^{1 - E'^){a - (3){r'± - I3)^^JT~~5. The moduli (variables) of the 
hypergeometric function of Appell are given by Eq.(36). 

Likewise we integrate the rest of terms in equation (61) and obtain equation 
(58) for the orbital period of a non-circular equatorial orbit. 



A. 0.2 Calculation of the angular integrals 

In this subsection we present some details of the computation that yields the 

polar orbital period (33). 

The polar orbital period is given by the following definite integral 



cP 



J rp 



dr{r'^ + a^) 

a7b~' 



d^ cos^ 



(70) 



VQ Jo \/i - k'2 sin^ * 

where Pr = E{r'^ + a^). The radial integral is calculated in a similar manner as 
in the previous appendix. Below we give the details of the angular integration. 
Using the latitude variable ^ the angular integral becomes 



-a'E 



ci*cos2(*) 



/"^ d^cos 



sin2(*) 



-a'E 



d*(l - sin2(*)) 



1-k'2 sin2(*) 



(71) 



Now the first angular integral using the variable x = sin \E' becomes 



-a^E 



d^ 



^ -'0 ^l-K'2 sin2(*) 



-a^E 

VQ 



dx 



V'(l-a;2)(l-Ac'2a;2) 



(72) 



Now using a new variable z defined by -^^r^x^ = z our angular integral 
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becomes 



-a2£;sin^^ /lllS.o ,o.o\^ -a^i; . ^/1113.o ,0.0 
i*! -,sin (p, «; sin 2 = — sin <^ i*! I sm y>, «; sin 



2 V2'2'2'2' ^' J VQ V2'2'2'2 

This result is used to get (31). Likewise the second angular integral is calculated 
elegantly in terms of generalized hypergeometric functions of Appell 

f"^ d^'(-sin2*) 



sin2(*) 

sin^ if f 3 1 1 5 .2 ,2 - 2^2 



„ Fi { 7:, 7:, 77, 7:, sin ifi,K sin oj - 
2 V2 2'2'2' ^' 3 

sint/jFi ( i i |,sin2(/3,K'2sin^</9] -siny-Fi ( ^ i i | , sin V, k'^ sin^ </? 



1 

772 '''M2'2' 2' 2'""" J '''"V2'2'2'2 

(74) 

where in the calculation the following values of the Gamma function have been 
used: T (i) = V^,T {^) = 0F/2,r(|) = ^,r(l) = 1. Use of radial inte- 
grations such as (65), (67), (68), (69) and angular integrations (73), (74) yields 
(33). 

B System of differential equations of Lauricella's 
multivariable hypergeometric function 

In this subsection we give some properties of the Fjj function introduced in (43). 
The fourth Lauricella function of m-variables is given by [37] 

77 /„ /J „^ _ \ ^ (<^)niH hwrn 

FD{a,p,j,z) - ^ z, ■■■z^ 

^ E(Qi,m H nm){f3i,n{) ■ ■ ■ {I3^,nm)zi^ ■ ■ ■ z'!^ 

^ „ {l,n\A \- nm)n\\ ■ ■ ■ UrJ 

where 

Z — (-21, ■ ■ ■ , ^m) 

/3 = (/?i,---,/3J (75) 
The Pochhammer symbol {a)m = {ct, m) is defined by 

r(a + m) J 1, if m = 0, 



r(a) I a(a+ l)---(a + rn- 1), if m = 1,2,3, • 
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It satisfies the following system of differential equations 

- Pi^l-S Pi^i-^a 



- PiZn^^ - af^Fo = (i = 1,2, • • • ,n) 



The series admits the following integral representation 



(76) 



Fn{a,^,r,z) = TTT^ r ['e-Hl-ty—\l-zrt)-('^---{l-Zmt)-''-dt 

r(a)r(7 - a) Jo 

(77) 

which is valid for Rc(a) > 0,Re(7 — a) > 0. It converges absolutely inside the 
m-dimensional cuboid: 

\zj\<l, (i = l,---,m) (78) 

C Exact orbital solution for non-spherical non- 
polar Kerr geodesies 

C.l Contribution to the change in azimuthal angle from 
radial integration for timelike orbits with L 7^ 

Using similar techniques as in section 2 we can calculate the contribution to 
the change of azimuthal angle from the radial integration for timelike non- 
spherical orbits with L ^ 0. The corresponding radial integral is Acj/" := 

2 J^^ aVr — where the quartic polynomial R is defined in (5) for 

A= 0. The above radial integral can be calculated exactly in terms of Appell's 
first generalized hypergcomctric function Fi of two variables 

- ^[ H+ ^i(^2'^'2'^'''+'^ j r(2) 



H- 2' ' -'^y r(2) 



H. ^^V2 2' ' r(i) 



(79) 
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where now the coefficients A^^ are given by the expressions 



iJVP 



r_ — r+ 



(80) 



C.1.1 Orbital solution for timelike non-spherical orbits with L ^ 

In this subsection and assuming a vanishing cosmological constant, we shall 
derive the exact solution for a non-spherical bound timelike orbit with non- 
vanishing parameters L and Q, i.e a non-spherical non-polar and non-equatorial 
orbit. The relevant differential equation is 



de 
7i 



(81) 



where now the polynomials R{r) and 0(6) are defined in (5) for A= 0. Using 
similar steps as in section 2.0.2 we obtain the exact orbital solution in terms of 
the square of Jacobi's sinus amplitudinous function 



BH 



a— 7 



1 _ £iz:^sn2 ( ^''^^y/i^-'rW-S) r de_ ^i 



(82) 



while the angular integration satisfies the equation 

^-a2(-l + i;2) + L2 + Q + ^a^{-\ + E^^Y - 2a2(-l + E-^){I? - Q) + {p VQf 



2 X 



4 1 1 1 ^ 5 



2^2 



where the modulus of Gaufi's hypergeometric function is given in terms of the 
roots of the quartic polynomial i?(r) by the expression 



2 _ a-/3(5-7 



a — 7 J — /? 

Solving for the angular variable ^ we obtain the exact expression 



(84) 



4 



sin^ * 



= sn 



62 - 63 



2\/ei - 63 



1 



^JY^rEi^f^r^4J^2 V2'2 



(85) 

where the three roots e, and the modulus k^" are given by equations (82), (83) 
in [11]. 

For i = 0, i.e. the case of polar orbits equation (85) reduces to equation 
(29). 



(83) 
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D Orbital period for equatorial non-circular Kerr 
geodesies in the presence of the cosmological 
constant 

The relevant differential equation is 

c'^' - ' "'("'+^') [ir'+a')E-aL] 1 

'dTr'Tn a;^ 71""^" ^ ^ ^''^ 

where the sextic polynomial is defined in (5) with Q = 0. Then the orbital 
period P^is given by the following integral 

''-''L '-';^°=<''^-^'*"^/, 7%- — -^r^ ^* 

(87) 

The first integral is given by the following exact expression 



oE^iaE-L) 

V w 



2' 2' 2' 2' ' ' ; r(i) 



..,3111 o 2 ■2\ r(i)r(i) 

+-^-(2'2'2'2'2'«'^'^j^|^ 



(88) 



where H,ijO and the arguments of the generalized hypergeometric function of 
Lauricella are given by equations (44), (45), (46) respectively. Now we calculate 
exactly the second term in the integration. First we have the integral 



r--aLHV±«2* (89) 

Jrp \/RAr 



Using partial fractions 



aLS2(r2+o2) A^' A^' A^' A^' 

^ ' - - + + — — + — (90) 



Ar ^ ~ ''a ^ ~ '^a r — r+ r — r 

the coefficients A*', i = 1, • • • , 4 are given by 



A'' 

A^' 



(ai - 0i){ai - 7i)(q!i - ^i)A 



(ai-(5i)(-/3i+(5i)(-7i+(5i)A 
^3, ^ 3(a3LS^ + «L/3^S^) 

(-ai+/3i)(-/3i+7i)(/3i-5i)A 
3(a3LS2 + aL7?H2) 



4/ 



(ai -7i)(/3i -7i)(7i -<5i)A 



(91) 
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Thus the integral I2 in closed analytic form will be given by generalized Lauri- 
cella hypergeometric functions Fd of four variables similar to those in eq.(50) 
and (53). 

In a similar fashion, the partial fraction expansion coefficients A^" for the 
third integral 

''^ ■E?a\r^ +<,^)E 



h :=2 

are given by the expressions 

A^" - 
A^" 

A^" : 

A^" 



-dr 



(ai - - 7i)(ai - ^i)-^ 

3{a^EE^ + a^ESlE^ 
{ai - + (5i)(-7i + ^1)^ 

3(a^^S^ + a'^E/SlE'' 
(-ai + /3i)(-/3i+7i)(/9i-5i)A 
3(a^^S^ + a^E-flE^ 
(ai-7i)(/3i-7i)(7i-'5i)A 



The fourth relevant integral is 

rrA E^rHr^+a^)E 
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VRAr 



dr 
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dr 



The first term in 1^ is calculated to be 



Jrp A y/R 



2j^E^E 



'U! 



-Fd 



1111 
2' 2' 2' 2 



.3111 2\2 2 

+ '^Pd (-,-,-, -,2,k , a ,/x 
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r(i/2)^ 
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(92) 



(93) 



(94) 



(95) 
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The partial fraction coefficients for the second term in are 

(?{a\ - /3i)(q!i - 7i)(ai - ^\)^ 
9{a^c^EE^ - 2^GMbh^iS^ + +c^EdlE^) 

c2(ai-^i)(-/3i+<5i)(-7i+<5i)A2 
9(a^c^j;S^ - 2EGMbiiI3iE^ + +c^^;/3iS^) 

c2(-ai+/?i)(-/3i+7i)(/3i-(5i)A2 

9(a^c^£;s2 - 2^GMbh7iS^ + c^^t^S^) 
c2(ai-7i)(/?i-7i)(7i-'5i)A2 



til// 



(96) 



E Periapsis advance for the orbits of S-stars 

In this appendix in order to gain an appreciation of the magnitude of the pe- 
riapsis advance for the S-stars wc perform a precise calculation of the effect a) 
using the exact formula for equatorial non-circular orbits produced in [12] and 
b) using the exact, closed form formula, equation (29), for non spherical polar 
orbits derived in this work. In the former case, wc assume that the angular 
momentum axis of the orbit is co-aligned with the spin axis of the black hole. 
In this case the effect of the spin of the black hole on the periastron precession 
is maximized. 

For the convenience of the reader we exhibit below the exact formula derived 
in [12], which is used in the first part of the current calculation 

A<^g™-27r = 



+ 




The derivation of the above formula as well as the definition of the coefficients 
A± and the cubic polynomial whose roots are u- can be found in [12], pages 
4412-4413. 

Our free parameters are the Kerr parameter a and the constants of integra- 
tion L, E. For a fixed initial choice of values for L, E and for different values of 
the Kerr parameter which include those supported from recent observations we 
calculated precisely the resulting periapsis advance, using the exact equations 
(97). The results are displayed in tables 14,15, 16. For comparison we note that 
for Mercury the perihelion advance for a particular choice of initial conditions 
and neglecting the rotation of the Sun (see table 1 in [6]), is: 0.103517 revolution • 
Thus we conclude that the theoretically calculated values for the periastron ad- 
vance for the S-stars, using (97), displayed in tables 14,15, are 521-7003 times 
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bigger than that of Mercury after one complete revolution. From the results of 
tables 14,15 we can see the effect of the rotation of the black hole. The effect 
is small and is of the order of a few arcseconds per revolution, maximized for 
the stars S2 and S14. From tables 14 and 15 we see that S2 and S14 come 
closer to the black hole than any other S-star at periastron passage. The total 
relativistic periapsis advance decreases by increasing a for fixed values of L, E. 
We repeated the analysis for different values of L, E. The results are presented 
in tables 17 and 18. Our results for the periapsis and apoapsis distances are in 
agreement with the data in [18]. 

Let us compare now the results for the periapsis advance in Tables 14,15 
with the predictions using the exact formula, equation (29). Indeed for the 
star S2, we calculated: A*<^™ ^ 683.916^,683.894^ respectively for a = 

7 rev ' rev ^ ^ 

0.52, 0.9939 for the values of Q, E in Tables 1,2. For S14 we derived the periapsis 
advance A^^itr = 731.855^^731.83^ respectively for a = 0.52,0.9939 for 
the values of Q, -B in Tables 1,2. This result is expected since for polar orbits the 
contribution of the spin of the black hole is minimal on the effect of periapsis 
advance. Our results of the calculation of periapsis advance using the closed 
form formula (29), for the S-stars are displayed in Table 19. The predictions for 
the periapsis and apoapsis distances are the same as those in Tables 14,15. 

On the other hand, let us compare the results of tables for the periastron 
advance in Tables 17 and 18 with those obtained using (29). For S2 we derived 
^^GTR ^ 738 445 arcs 738 arcs f^j. ^ ^ q 52 0.9939 rcspcctivcly and for the 

rev ' rev ' ^ 

choice of values for the parameters Q, E as those displayed in table 5. For S14 
we derived A^^tr ^ 926.551^,926.511^ for a = 0.52,0.9939 respectively 
and for the choice of invariant parameters Q, E for S14 as those exhibited in 
table 5. " 

The astrometric capability of GRAVITY will also allow deriving very accu- 
rate orbits for the S-stars. As we saw in the main text the orbital periods in 
this case are in the range of 14.8yr-94yr, and thus more precise astrometric data 
from GRAVITY will allow the probing of general relativistic effects on a longer 
(decade) time scale. 



^^For non-spherical non-polar orbits with orbital inclination 0° < i < 90° one expects that 
the resulting periapsis advance has a value between the ones calculated using formula (29) and 
(97) (assuming the 3 orbital configurations, polar, non-polax and equatorial have the same 
eccentricity and semi-major axis). A detailed investigation of this matter is however beyond 
the scope of this paper and will be analysed elsewhere. 
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Star 



E 



2-K 



Periapsis 



3.17 X lO^cm 

1.82 X IQi'^ cm 
2.87 X 10^^ cm 
3.35 X 10^5 cm 
1.58 X 10^6 cm 
1.64 X 10^5 cm 



Apoapsis 



6.69 X l(F^cm 
2.74 X 10^6 cm 
7.58 X 10^6 cm 
6.49 X 10^^ cm 
3.65 X 10^^ cm 
5.22 X 10^6 cm 



SI 

S2 
S8 
S12 
813 
S14 



267.778473 

75.4539876 
96.1191430 
103.1599 
192.08634 
72.9456205 



0.999993919 

0.999979485 
0.999992385 
0.999991213 
0.99998856 
0.999988863 



54.0874 
677.643 
418.185 
363.183 
105.016 
724.911 



revolution 

arcs 
revolution 

arcs 
revolution 

arcs 

revolution 

arcs 

revolution 

arcs 
revolution 



Table 14: Periastron precession for the S-stars in the central arcsecond of the 
galactic centre. The; Kerr parameter of the galactic black hole is acaiactic = 
Q 52^i^ii. The values of the parameter L are in units of We assume a 

central mass Mbh = 4.06 x 10^ Mq. 



Star 



E 



2lT 



Periapsis 



3.17 X 10^** cm 

1.82 X lO^^ cm 
2.87 X 10^5 cm 
3.35 X 10^^ cm 
1.58 X 10^^ cm 
1.64 X 10^5 cm 



Apoapsis 



6.69 X 10^** cm 
2.74 X 10^^ cm 
7.58 X 10^6 cm 
6.49 X lO^'' cm 
3.65 X 10^^ cm 
5.22 X 10^6 cm 



SI 

S2 
S8 
S12 
S13 
S14 



267.778473 
75.4539876 
96.1191430 
103.1599 
192.08634 
72.9456205 



0.999993919 

0.999979485 
0.999992385 
0.999991213 
0.99998856 
0.999988863 



53.9597 
671.947 
415.429 
360.953 
104.670 
718.607^ 



revolution 
arcs 

revolution 
arcs 

revolution 
arcs 

revolution 

arcs 

revolution 

arcs 
revolution 



Table 15: Periastron precession for the S-stars in the central arcsecond of the 
galactic centre. The Kerr parameter of the galactic black hole is acaiactic = 
0.9939 ^^^^ . The values of the parameter L are in units of We assume 

a central mass Mbh = 4.06 x lO^M©. 



star 


case l:(a = 0.9953) A(/)^^'"- - 2tt 


case 2: (a = 0.99979) A^e^'^ - 27r 


SI 

S2 
S8 
S12 
S13 
S14 


53.9593 

revolution 

671-930 ^.^^ 
415.420 T^ 

revolution 

360.9469 
104.669 j,,^ 
718.589 

revolution 


53.9581 "'^^ 

revolution 

671.876 Tt- 

revolution 

415.394 Tf- 

revolution 

360.926 

revolution 

104.666 Tt- 

revolution 

718.529 

revnlutinu 



Table 16: Periastron precession for the S-stars in the central arcsecond of the 
galactic centre, using the exact, closed-form formula, equation (97), for two high 
values for the spin of the black hole. The values for the invariant parameters 
L,E are the same as those in tables: 14,15. We assume a central mass Mbh = 
4.06 X lO^Mfl. 
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star 


L 


E 


A(/)^'^^ - 27r 


Periapsis 


Apoapsis 


S2 
S14 


72.6190898 
64.8441485 


0.999979145 
0.999987653 


■ ^ revolution 

916.661 


1.68 X 10^-^ cm 
1.29 X 10^5 cm 


2.71 X 10^" cm 
4.73 X 10^^ cm 



Table 17: Periastron precession for the S-stars in the central arcsecond of the 
galactic centre for different values of the parameters L, E than those of table 
14. The Kerr parameter of the galactic black hole is acaiactic = 0.52^^2. We 
assume a central mass Mbh = 4.06 x IQ^Mq. 



Star 


L 


E 


A^^'' " - 27r 


Periapsis 


Apoapsis 


S2 
S14 


72.6190898 
64.8441485 


0.999979145 
0.999987653 


725.018^;,^ 
907.686 


1.68 x lO'*' cm 
1.29 X 10^^ cm 


2.71 X 10'** cm 
4.73 X lO^'' cm 



Table 18: Periastron precession for the S-stars in the central arcsecond of the 
galactic centre for different values of the parameters L, E than those of table 
14. The Kerr parameter of the galactic black hole is acaiactic = 0.9939^^^. 
We assume a central mass Mbh = 4.06 x IO^Mq. 



Star 


case l:(a = 0.52) A*"'^'^ 


case 2: (a = 0.9939) A*^'^'^ 


SI 
S2 
S8 
S12 
S13 
S14 


54.2277 T^ 

revolution 

683.916 

revolution 

421.218 T!- 

revolution 

365.636 T^ 

revolution 

105.396 T^ 

revolution 

731.855 Tt- 


54.2275 T^ 

revolution 

683.894 ^f'^". 

revolution 

421.209 Tf- 

revolution 

365.63^,^,11^ 
105.395^,^11^ 
731.83 T^ 

revolution 



Table 19: Periastron precession for the S-stars in the central arcsecond of the 
galactic centre, using the exact, closed-form formula, equation (29), for two 
different values for the spin of the black hole. The values for the invariant 
parameters Q, E are the same as those in tables: 1,2. We assume a central mass 
Mbh = 4.06 x IO^^Mq. 
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